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ABSTRACT 

We provide a comprehensive description and offer an explanation for the sizes of the faintest known 
galaxies in the universe, the dwarf spheroidal (dSph) satellites of the Milky Way and Andromeda. 
After compiling a consistent data set of half-light radii (^1/2) and luminosities, we describe the size- 
luminosity relation of dSphs by a log-normal distribution in r%/2 with a mean size that varies as 
a function of luminosity. Accounting for modest number statistics, measurement uncertainties and 
surface brightness limitations, we find that the size-luminosity relations of the Milky Way and An- 
dromeda dSph populations are statistically indistinguishable, and also very similar: their mean sizes 
at a given stellar luminosity differ by no more than 30%. In addition, we find that the mean size, slope 
and scatter of this log-normal size description of Local Group dSphs matches onto the relation of more 
massive low-concentration galaxies. This suggests that the stellar sizes of dSphs are ultimately related 
to their overall initial baryonic angular momentum. To test this hypothesis we perform a series of 
high resolution N-body simulations that we couple with a semi-analytic model of galaxy formation. 
These predict the same mean size and slope as observed in dSph satellites. At the same time, these 
models predict that the size-luminosity distributions for satellite galaxies around similar host-halos 
must be similar providing a natural explanation as to why the size distributions of Milky Way and 
Andromeda satellites are similar. Although strong rotation is currently not observed in dSphs, this 
may well be consistent with our angular-momentum-based explanation for their sizes if the disks of 
these galaxies have become sufficiently stirred through tidal interaction. 

Subject headings: galaxies: dwarf, Local Group, galaxies: structure 
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1. INTRODUCTION 



Q Galaxies with stellar masses of 10 8 — 10 11 ' 5 M have a 
well defined relation between their stellar mass and the 
characteristic radius of their stellar body. This relation 
is one aspect of well-established global galaxy-parameter 
relations, such as the fundamental plane or the Tully- 
Fisher relation (IKormendvl Il977t iTullv fe Fisher! [l977t 
iDiorgovski fc Davislll987HDressler et allll987l ). For disk 
galaxies, the size of the stellar body is generally related 
to the angular momentum created by torques produced 
during the hierarchical formation of a galaxy and the 
infall of material towards the center of the gravitational 
potential well. In-falling gas dissipates and cools, but can 
only collapse until an gular momentum prevents further 
central concentration jFall fc Efstath iou I98Cfa iMo et al.l 
119981: iShen et al.ll2003D . This implies a size scaling with 
stellar mass, with the latter relating to the halo size, the 
baryonic fraction and the halo spin parameter. 

This simple angular momentum description is mod- 
ified by various mechanisms of galaxy formation and 
evolution: major merging for massive early-type galax- 
ies, questions relating to angular momentum conserva- 
tion, bar instabiliti es (which can redistribute internal an- 
gular momentum: [ Atha nassoula 20031: IKormendvl 1 19891 : 
Ivan den Boschlll998l ). as well as star formation efficiency. 
Despite these complicating effects, however, a simple 
angular-momentum-based formulation has been shown 
to provide a plausible explanation for the sizes of most 
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galaxi es at both low (e.g-.IShen et aLll2003[ ) and high red- 
shifts (|Somerville et al.l 120081) . At the low mass end of 
the galaxy scale, however, the size distribution of dwarf 
spheroidal galaxies (dSphs) which contain at most a few 
million stars is poorly determined empirically and with- 
out a cogent explanation. 

The sizes of Local G r oup d Sphs came into focus when 
IMcConnachie fc Irwinl (|2006D pointed out that bright 
dwarf galaxies of M31 were larger by a factor of 2-3 com- 
pared to their Milky Way counterparts of the same lu- 
minosity. As we will show later on in the paper, such an 
apparent size difference is difficult to reproduce in a cos- 
mological framework if the stellar sizes reflect halo sizes 
or halo angular momenta: N-body simulations fail to ex- 
plain a size discrepancy in sub-halo populations through 
differen ces in either host halo co llapse-time or host halo 
masses. IPenarrubi a et all (|2008l ) investigated whether a 
varied amount of tidal stripping between the Andromeda 
and Milky Way satellites could produce the observed size 
difference. They found that tidal effects could cause the 
reported difference in the characteristic sizes, but that 
the tides would also simultaneously lower their surface 
brightness. On the contrary, M31 dwarf galaxies appear 
to be both larger and lower surface brightness than Milky 
Way dwarf galaxies, thereby challenging this scenario. 
One caveat of these simulation s is that they do not in - 
clude a baryonic disk. Later on, IPenarrubia et all ([20101 ) 
have shown that including such a disk can have a signif- 
icant impact on the evolution of satellite dwarf galaxies. 
As a consequence, a difference in the stellar mass of the 
Milky Way and M31 disks could be driving some of the 
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apparent differences between their satellite populations. 
This analysis does not, however, explicitly study the pos- 
sibility of induced changes in size of the simulated dwarf 
galaxies. 

There are several practical difficulties in determining 
the sizes of faint, low surface brightness dwarf galaxies 
(My > —16). These galaxies can only be studied in 
detail within the Local Group, where they can be re- 
solved into their constituent stars: a single dwarf galaxy 
can contain as little as a few tens of stars down to cur- 
rent survey limits and determining their half-light ra- 
dius, T^^j^equires a suitable mathematical approach 
(e.g. IMartin et aD2008t iBrasseur et alJl2011al ). In addi- 
tion, survey detection limits and, therefore, the incom- 
pleteness in dwarf galaxy searches, play an important 
role in how we observe and build up their size-luminosity 
relation. If detection limits are not taken into account, 
one will invariably find that fainter galaxies tend to be 
smaller if large, faint galaxies (which have a lower surface 
brightness) fall belo w the detectability threshold. 

Since the work of iMcConnachie fc Irwin! ([2006) , the 
number of known Local Group satellites has strongly in- 
creased, (e.g. the number of known Andromeda dSphs 
has more than tripled). Additionally, since this time, 
satellites have been discovered more than a factor of 10 
fainter in luminosity, and we have increased our under- 
standing of the surface brightness limits of both Milky 
Way and Andromeda surveys. 

In this paper we set out to statistically quantify for the 
first time the global size-luminosity distribution of dSph 
satellites around the Milky Way and Andromeda, using 
a mathematical description which accounts for the detec- 
tion limits and low-number statistics. We then proceed 
to address the following questions: Are there differences 
in the mean sizes of Milky Way and Andromeda satellite, 
as previously pointed out? How does the size-luminosity 
relation at the extreme low-mass end of galaxies fit into 
a cosmological context? Does an interpretation of angu- 
lar momentum setting the characteristic size of galaxies 
apply to the very low mass end, and does this fit with 
current observations? 

In section [2j we compile a state-of-the-art data set 
of the sizes and magnitudes of low luminosity galaxies 
around the Milky Way and Andromeda. We then develop 
and apply a technique to describe the size-magnitude 
relation for galaxies around Andromeda and the Milky 
Way in section [3] We explore any possible size differ- 
ences between the sub-populations before comparing the 
r 1 / 2 — My relation found here (mean size, slope, scat- 
ter) to that of more massive galaxies (jShen et al.l 120031 : 
section [4]) . We put this in the context of cosmological 
expectations in section [3] and finally conclude on the im- 
plications for explaining the sizes of the smallest galaxies 
in section [6l 

2. THE SAMPLE 

In this study we are interested in only the faintest 
satellite galaxies and, as such, the Milky Way and An- 
dromeda satellite systems are the only two laboratories 
available for studying in detail the size-luminosity rela- 
tion of dSphs. Our sample therefore consists of all de- 
tected satellites of the Milky Way and Andromeda within 
a stellar mass range equivalent to — 14 < My < —6. We 



note that throughout this paper we use absolute magni- 
tude and stellar mass almost interchangeably, since the 
stellar mass-to-light ratios of these dwarf galaxies have 
little scatter (they are mainly old, metal-poor stellar pop- 
ulations) and vary smoothly with stellar mass. 

As much as we can, we strive to construct a sample 
of homogeneously measured sizes and luminosities of all 
Milky Way and Andromeda dSphs which we correct for 
the effects of surface brightness limitations. 

2.1. The Sizes and Luminosities of Milky Way and 
Andromeda Satellite Dwarf Galaxies 

One of the difficulties of putting together a homo- 
geneous sample of Local Group dwarf galaxy proper- 
ties comes from the strongly heterogeneous values that 
can be found in the literature, measured using differ- 
ent data sets, and different techniques. We have there- 
fore tried as much as possible to restrict ourselves to 
a small set of comprehensive st udies. Given these con- 
sidera tions , we re l y nam ely on llrwin fe~H atzidi mitrioul 
(| 19951) and lMateol (I1998D for the b right Milky Way dwarf 
galaxies, and IMartin et all (|2008f) for the recent Milky 
Way dwarf galaxy discov eries. For bright And r omed a 
dwarf galaxies, we rely on IMcConnachie fc Irwinl (2006), 
updated with measurements from dwarf galaxy discov- 
ery papers fro m the Pan- Andromeda Archaeological Sur- 
vey (PAndAS: IMartin et al.|[200l [20091: llbata et al.l l2007l: 
IMcConnachie et al.ll2008t IRichardson et al.ll201lD . a pho- 
tometric survey of the M31/M33 group conducted with 
the MegaCam wid e-field camera on the Can ada-France- 
Hawaii Telescope ([McConnachie et al.l [20091) . The PAn- 
dAS data themselves have been partially superseded by 
deeper photometric follow-up data sets tha t we have ob- 
tained for a large fraction of these satellites ([Collins et al.l 
2010; Bra sseur et al.l [2011 al Tbh . Detailed references are 
listed in Tables ED and H 

For dwarf galaxies with My < —8, size and luminos- 
ity measurements are quite robust and depend little on 
the method used to obtain them. Dwarf galaxies with 
My > —8, on the other hand, usually contain only a few 
tens or hundred stars down to the observational limit 
from which one has to reliably measure a size and a lu- 
minosity. Some of the difficu lties inherent to thi s prob- 
lem are described in detail in IMartin et al.l ((20 08) which 
also presents a maximum likelihood algorithm to mea- 
sure the structural parameters of a dwarf galaxy. We 
rely as much as we can on publications that have used 
this technique to derive the size of the galaxies in our 
sample. For details on its ap plications, the reader is re- 
ferred to IMartin et~aTl (120081) for the study of faint Milky 
Way satellites and IBrasseur et al.l (|2011al ) for an appli- 
cation to faint Andromeda satellites. 

An additional difficulty comes from the use of different 
effective radii in the literature. We have homogenized 
the data set by using, whenever possible, half-light radii 
measured along the major axis of the systems, and with 
the assumption that the radial density (or light) profile 
of the galaxy is exponential. When only the exponential 
scale radius is available, it has been converted to ri/ 2 
by multiplying by 1.68. Values that correspond to the 
geometric mean of the minor and major axes half-light 
radii have been converted using the ellipticities given in 
the corresponding study. 

For estimating total stellar luminosity, IMartin et al.l 
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Fig. 1. — Comparing the observed size distribution, ri/2, of An- 
dromeda low luminosity dwarf satellites galaxies (blue triangles) 
with those in the Milky Way (orange circles). The shaded regions 
correspond to the approximate observation limits of the Milky Way 
(dark grey) and Andromeda (light grey). 



(2008) have shown that color- magnitude diagram shot- 
noise can become an issue but this is only the case for 
very faint systems, so not essential for Andromeda satel- 
lites, which all have My < — 6.0 

2.2. (In) completeness Limits of Dwarf Galaxy Surveys 

All galaxy surveys suffer from some level of incomplete- 
ness; in the case of Local Group galaxies this arises from 
areal coverage, volume probed and the presence of Milky 
Way foreground stars. For example, the Sloan Digital 
Sky Survey, SDSS, led to the discovery of many new 
dSphs, but, as the survey area samples only one quar- 
ter of the sky around the North Galactic Cap, many of 
the known Milky Way satellites are concentrated in this 
region of the sky and undiscovered dSphs are likely to lay 
in regions of the sky which have not been imaged with 
the depth of the SDSS. 

Volume incompleteness issues do not, however, affect 
our final result if we assume the size-luminosity relation is 
not radially dependent. For our analysis we do, however, 
need to properly account for the luminosity and surface 
brightness limits of dSph searches around the Milky Way 
and Andromeda galaxies. 

In the case of the Milky Way, iKoposov et all (pOOc? ) 
quantify the distance dependent surface brightness lim- 
its for discovering satellites in the SDSS. This broadly 
corresponds to a limit of ~ 30 mag/arcsec 2 . In the 
case of Andromeda, all the dwarf spheroidal galaxies 
lie at approximately the same distance from us, lead- 
ing to a fixed surface brightness limit for discovering 
dwarfs. PAndAS has completed observations out to 150 
kpc from the center of M31 to an average surface bright- 
ness of ~ 29.2 mag/arcsec 2 (the central surface bright- 
ness of And XIX, the dwarf galaxy with the lowest sur- 
face brightness found in PAndAS), shown in light grey in 



Figure [T] Additionally, those galaxies with total magni- 
tudes fainter that My > —6 would not be easily detected 
in PAndAS, since they would have only a few stars along 
there red giant branch, and would go undetected by the 
selection techniques currently employed to discover new 
dwarfs. 

The PAndAS collaboration is currently in the process 
of accurately determining the completeness limits of their 
dwarf galaxy searches but this study can afford to rely 
only on these current rough estimates. We will show 
later that shifting the assumed detection limits for both 
Milky Way and Andromeda surveys to more conservative 
values does not change our results significantly. 

3. DETERMINING THE STELLAR 
SIZE-LUMINOSITY DISTRIBUTION 

In this section, we outline our method for determining 
the size distribution of dSphs as a function of magnitude. 
Any practical approach must account for the uncertain- 
ties, modest number statistics, and the observational de- 
tection limits of each sample. We satisfy these require- 
ments by adopting a parametrized functional form for the 
distribution and use a maximum likelihood algorithm to 
identify sets of parameters that make the observed data 
likely 

For our analysis we adopt a log-normal distribution 
for the sizes of dSphs. This is found to approxi- 
mate the intrinsic size distribution of mo re massive low- 
concentration galaxies (Shen et al. 2003) and, as will be 
detailed in section 14.21 is also well motivated by theoret- 
ical considerations: for massive disk galaxies, the size of 
the stellar disk is generally related to the angular momen- 
tum, with the latter relating to the halo spin parameter. 
The log-normal distribution of A leads to a size distribu- 
tion that is approximately log-normal. 

We fit the Milky Way and Andromeda dSphs sepa- 
rately in [logri/2, My]-space (seen in Figure [T| , adopt- 
ing the functional form from iShen et "all 120031 : a nor- 
mal distribution in log 7^/2 that has a slope, S, with 
magnitude and a dispersion, o~i gri around this slope. 

With lg r being the mean of logr 1 / 2 at My = —6.0 
and AMy = My + 6.0 being the magnitude offset from 
My = —6.0, we parametrize the mean of the distribu- 
tion, Is 
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A critical part of our analysis is accounting for the sur- 
face brightness limitations of each distribution. These 
limits are parametrized at the faint end as My (log 7*1/2), 
and are represented by the shaded areas in Figure [TJ 
To properly account for these limitations, we integrate 
our likelihood function over the observable window in 
[logr!/2, My]-space in order to correctly recover the nor- 
malization constant, A, of our likelihood function. Thus, 
A is re-computed for each choice of lg r, o~\ g r and S in 
the maximum likelihood grid, and we fix the normaliza- 
tion constant such that: 
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where n is the total number of dwarfs in the observable 
window. 

We also want to account for the individual measure- 
ment errors of each dwarf. To do this, we parametrize 
the variance of our distribution, cr[? r , as the sum of the 

intrinsic variance of the log- normal distribution, of r , 
and the error on each individual logri/a measurement, 

^lg r, i • 
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which reduces to erji r in the limit of zero measurement 
errors. Therefore the likelihood of a set of n dSphs, with 
the log of their half-light radius Xi, the corresponding 
uncertainties o\ g r , i, magnitudes My,i such that AMy,i = 
M Vi + 6.0, can be expressed as: 
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In practice it is often more convenient to work with the 
logarithm of the likelihood function. Since the logarithm 
is a monotonically increasing function, the logarithm of a 
function achieves its maximum value at the same points 
as the function itself. Thus our likelihood function which 
we compute at each position within a fine grid in param- 
eter space becomes: 

In C ({x l , trig r<i , AMy,J|lg r,o\ g r ,S) 

n 

= n In A(IgT, <7 lg r , S) - | ln(27r) - ^ In ( J crf g ~ + ^ 
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From there, we determine the model that maximize the 
likelihood function by exploring a fixed grid of parame- 
ters lg r, S and ai g r ,i between 1.99 and 2.55, 0.01 and 
0.60, and -0.4 and 0.0, respectively, and with step sizes 
of 0.01 in all cases. 

4. THE DWARF GALAXY SIZE-LUMINOSITY 
RELATION OF MILKY WAY AND ANDROMEDA 
SATELLITES 

We show the results of our likelihood method in Figure 
[2] for Milky Way (orange) and Andromeda (blue) dSphs. 
For these results, we have restricted the maximum likeli- 
hood method to those satellites with My < — 6 to ensure 
consistency between the Milky Way and Andromeda. 
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Fig. 2. — Parameters of a log-normal size distribution for the 
Milky Way (orange) and Andromeda (blue) dSph populations, de- 
rived from our likelihood analysis. The distribution of dSphs is as- 
sumed to be log-normal in i"]^, with a mean size at My = —6, lg r, 
an intrinsic dispersion, ci g r , and a slope in magnitude, S. Upper 
Row: The two dimensional relative likelihoods of each pair of the 
three fitted parameters marginalized over the third model parame- 
ter. The contours correspond to the 1 and 2-a limits, derived from 
X 2 values with 2 degrees of freedom: —2 ln£ r = 2.30 and 6.17, 
respectively. The most likely values are shown by dots in yellow 
(Milky Way) and cyan (Andromeda). Although the maxima for 
the Milky Way and Andromeda occur at different points, there is 
significant overlap in the 1 and 2-sigma contours. Thus, our global 
analysis finds no significant difference between the Milky Way and 
Andromeda dSph populations. Middle Row: The relative like- 
lihood functions of lg r, S and cri g r , marginalized over the other 
two model parameters for Milky Way dSphs. The colored regions 
under the curves indicate the 1-ct confidence intervals derived from 
X 2 values with 1 degree of freedom: —2 ln£ r = 1. Our computed 
values for each parameter are printed at the top of each window. 
Lower Row: Same as the middle row but for Andromeda dSphs. 



4.1. Comparing the Milky Way and Andromeda 
satellite systems 

In the top panel of Figure [D we overlay the result- 
ing 1 and 2-cr contours of the Milky Way (orange) and 
Andromeda (blue) size-luminosity parameters. Yellow 
and cyan dots correspond to the parameter values of 
the preferred models, i.e. the set of model parameters 
that maximizes the likelihood function for the Milky 
Way and Andromeda, respectively. The 1 and 2-o lim- 
its are derived from \ 2 values with 2 degrees of freedom: 
—2 ln£ r = 2.30 and 6.17, respectively. One can see from 
the top panel of Figure [21 that although the preferred 
models for the Milky Way and Andromeda are different, 
there is significant overlap of the 1 and 2-sigma contours. 

The lower two panels of Figure [5] show the resulting 
relative likelihood, C r = C/C ma x, functions for lg r, S 
and dig r , marginalized over the other two model param- 
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eters. The colored regions under these curves indicate 
the 1-er confidence intervals derived from \ 2 values with 
1 degree of freedom: — 2 ln£ r = 1. The Milky Way 
and Andromeda values of Igr (2.38i°;i3 and 2 - 35 ±ai4> 
respectively) and cr lg r (0.20t£u2 dex and 0.24±g;g| dex, 
respectively) are consistent within 1-cr. Only the deter- 
mined S values deviate slightly from each other, but only 
corresponds to a 1.7— a deviation for this parameter, with 
S = -0.03 ± 0.03 for the Milky Way and S 
for the Andromeda dSph population. 

One caveat of our results is that our maximum likeli- 
hood method relies on our choice of surface brightness 
limits. As described in section [2~2l these limits are not 
precisely defined. To test how the choice of these limits 
affects our results, we have additionally run our maxi- 
mum likelihood technique using more conservative sur- 
face brightness limits. We re-ran our algorithm for ten 
different choices of the surface brightness limits for the 
Milky Way and Andromeda, in steps of 0.2 mag/arcsec 2 
towards brighter limits. In each case, we again find sig- 
nificant overlap in the likelihood contours for Milky Way 
and Andromeda dSphs. We therefore conclude that for 
reasonable choices in the surface-brightness limits of the 
observations, our results remain unchanged. 

In addition, this analysis is performed assuming that 
two of the largest Andromeda satellites, And XIX 
and XXVII are bound systems, while observations of 
their surroundings show th at they are likely to be 
embedded in stellar st reams (jMcConnachie et ail 120081 : 
Richardso n et al. l l2011l ). This might hint that they are 
strongly deformed, and maybe inflated, by their interac- 
tion with Andromeda. Removing them from our sample 
would further diminish the small difference we measure 
between the size-luminosity relations of Milky Way and 
Andromeda dwarf galaxy populations. 

To conclude, we find that when one accounts for com- 
pleteness limits and the modest number of objects in 
the sample, there is no significant difference between the 
global size-luminosity relation of the Milky Way and An- 
dromeda dSph populations. 

4.2. The size distribution of the combined sample 

In the previous section, we find that the size-luminosity 
relations of Milky Way and Andromeda dSphs are sta- 
tistically indistinguishable. Building on this result, we 
assume a common size-luminosity relation for all dSphs 
by combining the Milky Way and Andromeda dSphs, and 
re-run our maximum likelihood method. In the lower two 
panels of Figure [3] we show the resulting marginalized 
relative likelihood curves with the colored regions under 
these curves indicating the 1-er limits of each parame- 
ter. The favored values for the model parameters are 
Igr = 2.39±g;g|, S = Q.24±S-^ and cr lg r = 0.06ig;g| dex. 
In the top panel of Figure y] we plot the resulting 1 and 
2-cr contours for these parameters, and the light purple 
dots correspond to the most likely model. 

4.3. A comparison with the sizes of more luminous 
galaxies 

Leaving the boundaries of the Local Group, one can 
wonder how the size-luminosity relation of dSphs com- 
pares to that of more massive galaxies. From a sample 
of 140,000 galaxies in the SDSS. IShen etail (|200l study 
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Fig. 3. — Same as for Figure [2] but with our maximum likelihood 
method applied jointly to all Milky Way and Andromeda dSphs. 

the size distribution of galaxies and its dependence on 
luminosity. They use a Sersic index of n — 2.5 to sep- 
arate their sample into high-concentration, early-type 
(n > 2.5) and low-concentration, late-type (n < 2.5) 
galaxies. 

The characteristic sizes of massive disk galaxies are 
generally related to the angular momentum of the sys- 
tem which is built up through torques during the process 
of hierarchical formation of a galaxy. Gas infalling into 
the gravitational potential well dissipates and cools, but 
is prevented from furt her collapse by the angular mo- 
mentum of the system dFall fc Efstathiou|[r980HMo et all 
1998; Shc n et al.l 120031 ). Therefore the gas, and conse- 
quently the stars, will retain a size scaling with the initial 
angular momentum of the system. Treating dark halos 
as singular isothermal spher es, and neglectin g the grav- 
itation effects of the disks, IMo et aD (|1998T ) derive the 
stellar disk scale length, Rd, to be: 



Rd 



~7k ( — I ^ r 200/r(A,Cvir,TOd, jd), (6) 

where ma and jd are the fractions of mass and angular 
momentum of the halo that are in the disk, A is the 
halo spin parameter, r2oo is the radius within which the 
mean density is 200 times the critical density, and f r is 
a factor that depends on the halo profile and the action 
of the disk: 
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Here, c v i r is the halo concentration (c V i r = r2oo/r s , with 
r s being the scale radius of the halo). 

N-body simulations show that the distribution 
of the hal o spin parameter t o be approximately 
log-normal dWarren et all Il992t ICole k, Lacevj 119961: 
Lemson fc Kauffmannl Il999t lAntonuccio-Delogu et al.l 
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Fig. 4. — The solid purple line shows our independent maximum 
likelihood fit to all Milky Way and Andromeda dSphs , while the 
dashe d blue and dashed pink lines correspond to the IShen et al.l 
(2003) relations for late- type (low-concentration) and early- type 
(high-concentration) galaxies, respectively. The shaded regions in- 
dicate the intrinsic scatter derived for both our work and the Shen 
et al. relation for late-type galaxies. 



P (\)d\ = 



27TCTi n . 



■ exp 



In 2 (A/A) 



This means that under the assumption of a constant 
jd/md the log- normal distribution of A will lead to a size 
distribution that is approximately log-normal. 

Motivated by these theoretical considerations, Shen et 
al. adopt a log-normal size distribution at a given lumi- 
nosity, characterized by dispersion ai g r . For faint, low- 
concentration systems, their results give a size depen- 
dence on r-band magnitude of \ogr 1 / 2 oc — 0.104 M r and 
dispersion, a\ g r = 0.19 dex. At the bright end, late-type 
galaxies have logri/2 oc — 0.204 M r and a\ g r = 0.12 dex. 
For comparison, these authors find that bright (M r < 
— 18.5) early-type galaxies follow logri/ 2 oc — 0.26 M r 
(see their equations (14), (15) and (16) for the full be- 
havior of their functional fits)Q 

In Figure |4] we plot our derived relation for dwarf 
galaxies (solid pu rple line) , with the relations derived by 
IShen et all ()2003| ) for brighter low-concentration (dashed 
blue line) and high-concentration (dashed pink line) 
galaxies in the SDSS. The shaded blue region indicates 
the intrinsic scatter in the derived relation for low- 
concentration galaxies in both our, and the Shen et al. 
relation. Remarkably, the mean size, the slope, and the 
dispersion found in our completely independent fit to 
the Milky Wa y and Andromeda dSphs, match beauti- 
fully on to the IShen et al.l (|2003l) relation for more mas- 
sive low-concen t ration galaxies. It should be noted that 
the IShen et aLl (|2003| ) low-concentration sample actu- 
ally includes dwarf elliptical galaxies (dEs) which, like 
dSphs, are low-concentration systems (e.g. Figure 11 

1 It should be noted her e that we have not attempted to convert 
the r-band magnitudes of Shen et al. (2003]) into V-band magni- 
tudes. Doing so may shift the Shen et al. relation by a small 
amount along the M^-axis, but not so significantly that it would 
impact our conclusions. 



of IGraham fc Worlevl 12008). The size-luminosity rela- 
tion we measure in the Local Group for dSphs galaxies 
is therefore also a continuation of that of dEs, hinting at 
a connection between these two types of dwarf galaxies. 

5. SIMULATION PREDICTIONS FOR THE SIZES 
OF FAINT DWARF GALAXIES 

We now present theoretical predictions for the size- 
luminosity relation of dSphs in a A Cold Dark Matter 
(ACDM) universe. These predictions are obtained by 
combining merger trees extracted from very high reso- 
lution N-body simulations describing the hierarchical as- 
sembly of a Milky Way-like halo, with semi analytic mod- 
els (SAMs) for galaxy formation. We perform ACDM 
simulations and analyze them, including the definition 
of dark matter halos and the reconstruction of their de- 
ta iled merger tree, usin g the tools described in full detail 
in lMaccio etaLl (|2010l ). hereafter M10. 

5.1. N-body code and semi analytical model 

We perform a new series of N-body simulations car- 
ried out using PKDGRAV, a treecode written by Joachim 
Stadel and Thomas Quinn ([Stadell 12001). with cosmo- 
logical parameters A = 0.742, £l m = 258, h = 0.72, 
n = 0.963 and er 8 = 0.79 (|Komatsu et al J 120091) . We se- 
lect nine candidate halos from an existing low resolution 
dar k matter simulation (350 3 particles within 90Mpc, 
see lNeistein et aLll2010l ) and re-simulate them at higher 
resolution using the volume re-normalization technique 
(jKatz fc White! 119931) . The high resolution run is 12 3 
times better resolved than the low resolution one: the 
dark matter particle mass is = 3.65 x 10 5 M Q , where 
each dark matter particle has a spline gravitation (co- 
moving) softening of 476 pc. 

To test the importance of the host galaxies' halo 
mass and formation time in shaping the properties of 
dwarf satellites, we select three halos for each of the 
following dark matter masses : Mi = 0.5 x 10 12 M Q , 
M 2 = 1.0 x 10 12 M o , and M 3 = 5.0 x 10 12 Af Q . For 
each mass we then select three halos of varied concen- 
tration: one with an average value (determined using 
the concentration-mass relation from M accio et al.l 120081 
and iMunoz-Cuartas et al.l 1201 ID and two with concen- 
trations that are 1.5-er higher or lower, than this aver- 
age. Since con centration strongly c orrelates with forma- 
tion time (e.g. iWechsler et al.ll2006l ) we therefore sample 
halos with early, average and late formation times for 
each mass bin. The properties of the individual halos 
are listed in Table 

In order to predict the expected luminosities of satellite 
galaxies, we combine the result s of the iV-body simula- 
tions with the SAMs outlined in lKang et all (|2005l I200H 
and M10. The effects of the various physical processes in 
shaping the luminosity function of Milky Way satellites 
were studied in full detail in M10. We determined that 
suppression of gas infall by a photo-ionizing background, 
supernova feedback and tidal destruction are the most 
relevant processes responsible for the agreement between 
theoretical predictions and observational data. We char- 
acterize our best fit models as follows: (i) we regulate 
star formation efficiency in low-mass halos by shutting 
off gas cooling in structures with virial temperature be- 
low 10 4 K (due to the inefficiency of H 2 cooling); (ii) we 
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Fig. 5. — Left: Sub-halo dark matter mass, M acc , and concen- 
tration, c v i r , at the time of accretion for a range of host halo masses 
and formation times (see Table [3] for details on each host galaxy). 
The upper plot shows three halos with the same formation epoch 
but different masses, and the lower plot shows three halos of the 
same mass, but different formation epochs with no noticible dif- 
ference in the satellite properties. Neither change in the host-halo 
mass nor changes in the formation time result in a systematic dif- 
ference in the concentration of the sub-halos. Note that Galaxy 1 is 
the highest mass, and hence shows the largest number of sub-halos. 
Right: The solid purple line and the shaded region correspond to 
our derived relation for the combined sample of Milky Way and 
Andromeda dSphs. The green and yellow points correspond to the 
sub-halos of simulations G4 and G5. 

suppress hot gas accretion in low-mass halos according to 
photo-ionization background, re-ioniza tion and filtering 
mass arguments (jOkamoto et al.ll2008| ): (iii) stellar feed- 
back is modeled as a function of halo circular velocity. 
We keep the parameters fixed to the values providing the 
best fit to observations in M1C0. None of these aspects 
were 'tuned' to match the sizes of the satellite galaxies. 

We first study to what extent the mass, M acc , and 
concentration, c vrr , of the sub-halos at the time of ac- 
cretion change as a function of their host halo mass and 
formation time. In the upper left panel of Figure El we 
plot the c v j r — M acc distributions for the sub- halo popula- 
tions of three host-halos with the same formation epoch, 
but different masses (0.41 x 10 12 M o , 0.77 x 10 12 M Q , and 
4.44 x 10 12 M Q ) and find no discernible difference between 
the three distributions. Similarly, we find no difference 
between the c v i r — M acc distributions of the sub-halo pop- 
ulations of three host-halos of equal mass but different 
formation epochs (see the lower- left panel of Figure [5]) . 

The similarity of these c vu — M acc distributions predicts 
indistinguishable satellite size distributions (see the de- 
pendence of Cvir and Af acc on the size of a stellar body in 
equation [6]). Therefore, within a cosmological framework 
we expect the same size-mass or size-luminosity relation 
for satellite halos orbiting host halos of different collapse 
times or massefl 

5.2. Comparison With Observations 

In order to compare our simulated satellites with ob- 
servations of dSphs, we first need to compute the size 

2 A more detailed and critical discussion of how the variation in 
SAM parameters may affect the final luminosity function can be 
found in M10. 

3 It should be noted that our simulations do not self-consistently 
include the presence of baryons in the main galaxy (e.g. a stellar 
disk), or in the satellites (a stellar core). While the presence of a 
stellar disk in the halo of the host gala xy can increase the mas s-loss 
of satellites with cored matter profiles (Pcnarrubia ct al. 2010), the 
satellite's ce ntrally concentrated stars will m ake it more resilient to 
tidal forces ( Schewtschcnko & Maccio 2011) and likely counteract 
this effect. Overall, it appears unlikely that the presence of baryons 
will dramatically change the results of our N-body simulations, but 
this should nevertheless be checked when moving beyond the simple 
model presented here. 



of the stellar body and the luminosity of each sub-halo 
within our simulation. While luminosity is an output 
of the SAM, the size of the stellar body, Rd, can be 
computed directly from N-body simulations according to 
equation To evaluate Rj(\, c vir , m d , j d , r 2 oo(-^haio)) 
we use the values of M acc , A and c v j r at the time at which 
the satellite was accreted since in our SAM star forma- 
tion is strongly suppressed after accretion (e.g. after a 
halo enters the virial radius of a more m assive halo). Ad- 
ditionally, we follow IShen et~aT1 (f200l in assuming that 
the material in a galaxy has specific angular momentum 
similar to that of the halo (i.e. jd/rrid = 1). 

In the right panel of Figure [S] we show the satellite 
populations of two simulated galaxies (yellow triangles 
and green circles) overlaying our derived relation for Lo- 
cal Group dSphs (solid purple line, with derived intrinsic 
scatter shown in blue) . For this comparison we have cho- 
sen two halos which are closest to the Milky Way in their 
properties, but, as we have shown above, there is no dif- 
ference in the distributions of our simulated sub-halos so 
this choice is arbitrary. One can see that the simulated 
galaxies match the mean size and slope of our derived re- 
lation well, with only the dispersion being slightly larger, 
stemming from the scatter in A. 

6. DISCUSSION 

We set out to describe the size-luminosity relation of 
Milky Way and Andromeda dSphs by a log-normal func- 
tion in ?T/2, characterized by an intrinsic dispersion, 

a\g r , and mean, lg r, modified to account for a slope, 
S, in magnitude. Using a maximum likelihood method 
which accounts for modest number statistics, observa- 
tional uncertainties as well as surface brightness limita- 
tions, we determine the size-luminosity relation of the 
Milky Way and Andromeda dSphs both separately and 
as an entire population. Our empirical results can be 
summarized as follows: 

1. There is no statistically significant difference be- 
tween the global size-luminosity relation of Milky 
Way and Andromeda dSph satellites. Indeed, the 
mean sizes of galaxies at a given stellar mass agree 
within 30%. 

Milky Way: 

log n /2 = 2.38±g;il - 0.03 ± 0.03 x (My + 6) (9) 
Andromeda: 

logr 1/2 = 2.35±8;E - 0-09t° °ol * {My + 6) (10) 

2. We find that the size- luminosity relation obtained 
for dSphs is in very good agreement with the rela- 
tion measure d for more massive , low-concentration 
galaxies from lShen et al.l (|2003f ). 

In a cosmological framework the first of our empir- 
ical results is comforting as our numerical simulations 
confirm that variations in host-halo mass or collapse- 
time/concentration do not lead to size differences for 
their sub-halos. This does not mean, however, that 
there do not exist Andromeda dSphs that are larger than 
Milky Way dSphs at a give n luminosity. These exist 
(|McConnachie fc Irwinl 120061: iMcConnachie et al.l 12001: 
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iMartin et all 120091 iRichardson et alj|2011l) but they do 
not translate into a significant difference in the global 
properties of the dSph populations. Our analysis shows 
that sampling, combined to the larger number of An- 
dromeda satellites, are enough to explain the presence of 
larger satellites around Andromeda. One does not need 
to invoke different formation mechanisms, or evolution 
pathways, to explain apparent size differences. Hints of 
dynamical differences between the memb e rs of the two 
satell ite populations (|Collins et al.l 201 0; K alirai et al.l 
1.201 Of ) could, however, be an tell-tale sign that the story is 
more complicated than displayed by their size-luminosity 
relation alone. 

With our dSph size-luminosity relation matching so 
well to that derived independently for more massive late- 
type galaxies, our second result indicates that there may 
exist a common size determinant for galaxies across a 
wide range of masses. For disk galaxies, the size of the 
stellar body is generally related to the angular momen- 
tum created by torques produced during the hierarchical 
formation of that (satellite) halo, and the infall of ma- 
terial towards the center of the gravitational potential 
well. If a fraction of the initial angular momentum is 
transferred from the halo to the baryons, the baryons 
proceed to cool into a disk until they are prevented from 
further collapsing by the angular momentum of the ma- 
terial. This results in a relation between the final size of 
the galaxy and stellar mass, which is itself related to its 
dark matter halo mass and the initial angular momentum 
o f the baryon s. 

IShen et all (pOOl were able to explain the size- 
luminosity relation of massive late-type (n < 2.5) galax- 
ies by applying the angular momentum formulation of 
I Mo et aU (|1998| ). Dwarf spheroidal galaxies, like dwarf 
elliptical galaxies, sh ow Sersic indices of n ~ 1 (e.g. 
iDe Rijcke et al.l[2003l ). making them a faint extension of 
late-type systems by this criterion. When we apply this 
simple angular momentum formulation derived for more 
massive disk galaxies to the scale of dSphs, we are able 
to reproduce the observed size-luminosity distribution of 
dSph galaxies in both mean size, slope and scatter (see 
Figure E|) . This strongly suggests that angular momen- 
tum arguments and the cosmological framework play an 
important role in setting the sizes of dSphs. 

At face value, such an explanation for the sizes of 
dSphs implies that they should be, or have been, rotat- 
ing. While rotation is not manifest in the existing data 
for many galaxies, there are signs of sig nificant net an- 
gular momentum in emergent data sets. [Battagli a et al.l 
( 2008) has found a radial velocity gradient in the Sculp- 
tor dwarf galaxy, whic h they interpret as an indication of 
internal rotation, and iFabrizio et al.1 ()2011h find the ra- 
dial velocity distribution across the body of the Carina 
dSph to show a radial gradient that they take as possible 
evidence of rotation. The explanation for these radial 
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velocity gradients is still debated, though, as they could 
simply be the result of the dwarf galaxies' proper motions 
(iKaplinghat fc Strigari| [2008; W alker et al.ll2008t iStrigarl 
I2010f ). 

Even if more subtle signs of rotation could be discov- 
ered in the future, the question arises of how clear sig- 
natures of rotation should be. The flatness and coldness 
of a gas disk, and the resulting stellar disk, is given by 
the ratio of its rotation velocity to its velocity disper- 
sion (v ro t/<7). Clearly, for low mass galaxies, this ratio 
will approach unity if the rotation velocities approach the 
turbulent velocities of the ISM, ~ 5kms _1 . In addition, 
much of the initial rotation in dwarf galaxies may have 
been erased if they have been tidally stirred via their in- 
te ractions with their host galaxy. The recent simulations 
of lLokas et al.l (|2011l ) postulate that such stirring indeed 
could or should be effective in erasing the rotational sig- 
natures in dSphs while, at the same time, only mildly 
reducing the size of their stellar c omponent. S u ch a co n- 
clusion also follows the work of iMaver et al.l (|2001alfbT l 
who initially proposed that dwarf irregular galaxies are 
actually the progenitors of dSphs, transformed through 
tidal stirring. This has some observational evidence, as 
it explains the fact that dwarf irregulars are found on the 
outskirts of their host galaxies, while dSphs are mainly 
located further in. 

The scenario we propose for the sizes of low-mass 
satellite galaxies would therefore start out with a well 
founded theory based on angular momentum considera- 
tions: dSph galaxies form as disks with sizes set by the 
angular momentum, and initially have v rot /cr ~ 1. Tidal 
stirring would then randomize some fraction of their kine- 
matics, leaving their disks with little sign of rotation, 
which is consistent with the data. In this theory, more 
distant satellites, which feel weaker tides should, on av- 
erage, show more rotation, as seems to be the case with 
the Cetus dwarf galaxy, an isolated Local Group dSph 
that shows a rotation signature of 8kms _1 (|Lewis et al.l 
I2007T ) . consistent with this scenario. 
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TABLE 1 
Properties of Milky Way dSphs 



dSph 



M 



v 



r l/2 



(pc)l 



Reference 



Bootes I 
Bootes II 
Carina 

Canes Venatici I 

Canes Venatici II 

Coma Bernices 

Draco 

Fornax 

Hercules 

Leo I 

Leo II 

Leo IV 

Leo V 

Leo T 

Sculptor 

Segue 1 

Segue 2 

Sextans 

Ursa Major I 

Ursa Major II 

Ursa Minor 



-6.3 ±0.2 
-2.7 ±0.9 
-9.3 

-8-618:? 

-4.9 ±0.5 
-4.1 ±0.5 
-8.75 ±0.05 
-13.2 
-6.6 ±0.3 
-11.9 
-9.6 

-5.8 ±0.4 
-5.2 ±0.4 
-8.1 ±0.1 
-11.1 
_i =+0.6 

-2.5 ±0.3 
-9.5 

-5.5 ±0.3 
-4.2 ±0.5 
-8.8 



2A2+_f 
51 ± 17 
241 ± 23 
564 ± 36 
74+" 

77 ±10 
196 ± 12 
668 ± 34 
330±g2 
246 ± 19 
151 ± 17 
206 ± 36 
133 ± 31 
178 ± 39 
260 ± 39 
29tg 
34 ±5 
682 ± 117 

31811° 
140 ± 25 
280 ± 15 



Martin et al. (2008) 
Martin et al. (2008) 
Irwin & Hatzidimitriou 
Martin et al. (2008) 
Martin et al. (2008) 
Martin et al. (2008) 
Martin et al. (2008) 
Irwin & Hatzidimitriou 
Martin et al. (2008) 
Irwin & Hatzidimitriou 
Irwin & Hatzidimitriou 
de Jone et al. (2010) """ 
de Jong et al. (2010) 
Martin et"al""(2008) 
Irwin & Hatzidimitriou 
Martin et al. (2008) 
Belokurov et al. (2009) 
Irwin & Hatzidimitriou 
Martin et al. (2008) 
Martin et al. (2008) 
Irwin & Hatzidimitriou 



1995), Mateo (1998 




1995), Mateo (1998 



1995), Mateo (1998 



1995), Mateo (1998 



Where 



1/2 



is the half-light radius measured along the semi-major axis. 



TABLE 2 
Properties of Andromeda dSphs 



dSph 



M 



V 



r l/2(pc) 1 



Reference 



Andromeda I 
Andromeda II 
Andromeda III 
Andromeda V 
Andromeda VI 
Andromeda VII 
Andromeda IX 
Andromeda X 
Andromeda XI 
Andromeda XII 



Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 
Andromeda 



XIII 

XIV 

XV 

XVI 

XVII 

XVIII 

XIX 

XX 

XXI 

XXII 

XXIII 

XXIV 

XXV 

XXVI 

XXVII 



-11.8 ±0.1 
-12.6 ±0.2 
-10.2 ±0.3 
-9.6 ±0.3 
-11.5 ±0.2 
-13.3 ±0.3 

8 1+0-4 

-7.36 ±0.07 

-6.9183 

-6.418;^ 

-6 7+°- 4 
°-'-o.i 

-8.46 ±0.07 

-9.8 ± 0.4 

-9.2 ± 0.5 

-8.30 ±0.06 

-9.66 ±0.07 

-9.57 ±0.69 

-7.16 ±0.10 

-9.33 ±0.07 

-6.2 ±0.15 

-10.2 ±0.5 

-7.6 ±0.5 

-9.7 ± 0.5 

-7.1 ± 0.5 

-7.9 ±0.5 



U3S 



*!! 4 r 
1126 ± 17 
403 ± 17 
336 ± 17 
454 ± 17 
739 ± 34 

= = 9+22 
ooz -110 

235 ± 20 

i40 -20 

28911? 
203±f 4 
406 ± 27 
270 ± 30 
136 ± 14 
265 ± 19 
317 ±23 
1647 ± 252 
160 ± 30 
777 ± 53 
405 ± 27 
1316 ± 75 2 
349 ± 19 2 
661 ± 37 2 
222 ± 12 2 
361 ± 21 2 



McConnac hie fc Irwin (2006) 
McConnachie fc Irwin (2006) 
McConnachie fc Irwin (2006) 
McConnachie fc I rwin (2006) 
McConnachie fc I rwin (2006) 
McConnachie fc Irwin" (2006) 
Collins et al. (2010) 
J3r^s^ejrr_et_aL (2011a) 
Collins et al. (2010) 
Collins et al. (2010) 
Collins et al. (2010) 
Brasseur et al. (2011b) 
Ibata et al""(2007) , Letarte et al. (2009) 
Ibata et al. (2007) . Letarte et al. (21)59*) 

Brasseur et al. (2011b) 
Brasseur et al. (2011b) 
Brasseur et al. (2011b) 
Brasseur et al. (2011b) 
Brasseur et al. (2011b) 
Richardson et al. (2011) 
Richardson et al. (2011) 
Richardson et al. (2011) 
Richardson et al. (2011 ) 
Richardson et al. (2011 ) 



1 Where rj/ 2 is the half-light radius measured along the semi-major axis. 

2 The uncertain t ies on these values are based on the distance uncertainties listed in 
Richardson ct al. (20TTD. 



TABLE 3 
Simulated Galaxy Parameters 



Halo 


Maw f 1 D 12 Mr- 


1 N 4.' 1 

J ly particles 




nV^TTTl atl(~lTl ("imO 


Gl 


AAA 


12,157,818 


339.5 


average 


G2 


4.46 


12,103,298 


339.9 


high 


G3 


4.89 


13,270,207 


349.8 


low 


G4 


0.77 


1,889,570 


189.8 


average 


G5 


0.89 


2,184,048 


199.2 


high 


G6 


0.81 


1,987,729 


192.2 


low 


G7 


0.41 


1,043,557 


154.4 


average 


G8 


0.48 


1,221,726 


161.3 


high 


G9 


0.51 


1,298,083 


165.0 


low 



